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^ I Abstract. We study the effect of non-trivial sound speeds on local-type non-Gaussianity 

lO ' during multiple-field inflation. To this end, we consider a multiple-DBI model and use the 

6N formalism to track the super-horizon evolution of perturbations. By adopting a sum 
\^ [ separable Hubble parameter we derive analytic expressions for the relevant quantities in the 

^f^ ' two-field case, valid beyond slow variation. We find that non-trivial sound speeds can, in 

^^ ■ principle, curve the trajectory in such a way that significant local-type non-Gaussianity is 

produced. Deviations from slow variation, such as rapidly varying sound speeds, enhance 
this effect. To illustrate our results we consider two-field inflation in the tip regions of two 
warped throats and find large local-type non-Gaussianity produced towards the end of the 
/\ , inflationary process. 
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1 Introduction 



Whilst inflation provides a compelling mechanism for producing the initial conditions on 
which the hot big bang relies, a consistent model proves elusive (see [1-3] for reviews). This 
is due, in part, to the limited information available from the two-point statistics of primordial 
density perturbations. Interest has therefore focused on potential non-Gaussian signatures to 
discern between otherwise degenerate models [4], particularly since the sensitivity of forth- 
coming observations is set to improve by at least an order of magnitude [5]. As such, it is 
important to understand the correspondence between inflationary dynamics and the resultant 
form of non-Gaussianity (see [6, 7] for recent reviews). 

Inflation is capable of producing non-Gaussianity in a variety of ways. For example, 
the conversion between entropy and adiabatic modes during^ multiple-field inflation [15-17] 
allows the curvature perturbation C to evolve on super-horizon scales [18, 19]. Such non- 
linearities can, in certain circumstances, produce local-type non-Gaussianity during inflation 
[20-42]. Some simplifying assumptions are usually made to make such calculations analyti- 
cally tractable however, owing to the increased phase space. For example, Vernizzi & Wands 
[21] calculated the non-linearity parameter in canonical two-field models possessing a sum 
separable potential V{4>,x) = V^'^'{(j)) + V^^'{x)- Byrnes & Tasinato [29] extended this be- 
yond slow roll by exploiting the Hamilton-Jacobi formalism [43, 44] and a sum separable 
Hubble parameter H{(j),x) = H^'^'{(f)) + H^^'{x)- The broad conclusion is that a turn in 
the trajectory, caused by the potential, can produce significant local-type non-Gaussianity 
during infiation. Furthermore, the effect is enhanced by violations of slow roll. 

Single field models with non-standard kinetic terms, often motivated by string theory, 
provide an alternative source of non-Gaussianity (see [45, 46] and references therein). The 
models we are concerned with have a characteristic sound speed Cg, where c^ = 1 in the 
canonical case.^ Non-Gaussianity is caused by quantum field interactions on sub-horizon 
scales and is typically of equilateral form. In the Dirac-Born-Infeld (DEI) case [48, 49], 
for example, a probe D-brane moving along the radial direction of a warped throat causes 
inflation. Here the radial co-ordinate of the D-brane plays the role of the inflaton. In this 
case the sound speed acts as an effective speed limit, facilitating inflation on a steep potential 
and producing significant equilateral type non-Gaussianity. 

In more general models it is expected that both contributions to non-Gaussianity will be 
relevant. The combined effect of non-standard kinetic terms and multiple- field dynamics has 
been considered by [50-60] , with emphasis generally on the effect of multiple-field dynamics 
on equilateral type non-Gaussianity. Here however, we focus on the effect of non-canonical 
kinetic terms on local-type non-Gaussianity, produced by multiple dynamical fields (see [61- 
63] for alternative examples). Specifically, we consider whether non-trivial sound speeds, 
as opposed to the potential, can produce local-type non-Gaussianity during a turn in the 
trajectory. 



^Alternatively, non-linearities can develop through the conversion between entropy and adiabatic modes 
after inflation. For example, the curvaton mechanism [8-10] and modulated reheating [11, 12]. See [13, 14] 
for recent realisations of such scenarios in string theory. 

^Whilst commonly referred to as the sound speed, it should be noted that this is technically the phase 
speed of fluctuations. See [47] for further discussion on this point. 



It should be noted that, in the canonical case, the effect is enhanced by deviations 
from slow roll. It is reasonable to expect similar behaviour in the scenario considered here, 
corresponding to rapidly varying sound speeds. Such behaviour has been explored by [64-68] 
in the single field case, who find additional features in the equilateral limit of the bispectrum. 
In this paper however, we show that rapidly varying sound speeds in multiple-field models 
can lead to large non-Gaussianity of local form. 

In this work we choose a multiple-DBI model, akin to that of [58-60], as a concrete 
example of multi-component inflation with non-standard kinetic terms. We then employ the 
6N formalism to track the super-horizon evolution of perturbations using the field fluctuations 
at horizon exit and the subsequent background trajectory. The trajectory is a solution of 
the homogeneous equations of motion, which we write in Hamilton- Jacobi form and derive in 
an appendix. Whilst the dynamics are described more naturally in this form, it also allows 
the adoption of a sum separable Hubble parameter, as in [29], to treat the two-field case 
both analytically and beyond slow variation. With these expressions, we ascertain whether 
non-trivial sound speeds can produce a turn in the trajectory during which local-type non- 
Gaussianity is produced. To illustrate our results we consider infiation in the tip regions of 
two warped throats, where one expects the sound speeds to increase rapidly and violate slow 
variation. 

The outline of the paper is as follows. We begin in section 2 by introducing the multiple- 
DBI model, before introducing and reformulating the relevant quantities using the 6N formal- 
ism. By utilising a sum separable Hubble parameter we analyse these expressions analytically 
in section 3, before summarising the necessary conditions for large local-type non-Gaussianity. 
We illustrate these results in section 4 by studying inflation in the tip regions of two warped 
throats. Finally, we conclude in section 5. 

Throughout this paper we use the (—,+,+,+) metric signature and set Mp = c = 1, 
where Mp = /— — ^ is the reduced Planck mass. Capital latin indices label scalar fields and 
any summation is explicit. Greek indices label space-time co-ordinates whilst lower case latin 
indices label spatial co-ordinates only, where the Einstein summation convention is adopted. 
Finally, commas denote partial derivatives and over-dots represent derivatives with respect 
to time. 



2 Multiple-DBI inflation, non-Gaussianity and the 6N formalism 

We begin this section by introducing a multiple-DBI model as a specific realisation of multi- 
component inflation with non-standard kinetic terms. With this model in mind we introduce 
the relevant quantities, paying particular attention to non-Gaussianity. Finally, we recast 
these expressions with the 5N formalism to study their evolution on super-horizon scales, of 
which we are primarily interested. 



2.1 Multiple-DBI inflation 

We begin by considering the following action 



ll'dS:V^[R + 2'^P,-2V], (2.1) 



* 2 



where Pj is a function of the single scalar field (/>/ and kinetic function Xj = —2g^^4'i,fj,4>i,u7 
whereas the potential y is a function of the set of scalar fields (j) = {</>i, </>2, •••, (/'at}- For the 
spatially flat Friedmann- Robertson- Walker (FRW) metric the collection of scalar fields act as 
a perfect fluid with isotropic pressure P = ^j Pj—V and energy density p = ^j- pi+V, where 
Pi = 2XjPj^Xi — Pi- In appendix A we derive the corresponding homogeneous equations 
of motion by recasting the action (2.1) using the Arnowitt-Deser-Misner (ADM) formalism 
[69]. This ensures that the results are in Hamilton- Jacobi form [43, 44] in which the Hubble 
parameter is a function of the scalar fields and takes precedence over the potential. Not only 
is this more suited to the case of non-trivial sound speeds but it will also prove useful for 
studying violations of slow variation. The resultant equations of motion are 



Pl,Xj 



H,i, (2.2) 



3^' = E (^ - P^) + ^' (2-3) 



where H is the Hubble parameter and is, in general, a function of the set of scalar fields 
Furthermore, it is convenient to define individual sound speeds 



ain = EMl = ?L^ . (2.4) 

We now model infiation as driven by A'^ probe D3 branes traversing A^ distinct warped throats 
glued to a compact Calabi-Yau manifold in type HB string theory.^ This was first considered 
by [58-60] to study the effect of multiple sound horizons on equilateral type non-Gaussianity. 
The corresponding expression for Pi is 



^/ = y(7)(^l-Vl-2/^'^^/J- (2-5) 

Here f^' parameterises the warped brane tension of throat / and is a function of 0/ only. 



'^This model is then distinct from the first example of multiple- field DBf [70], in which a single brane 
descends a warped throat along radial and angular coordinates. 



On substitution into (2.2)-(2.4) we find the following equations of motion 

<i)i = -2c^'^Hj, (2.6) 



^'''-y-E^i^-j, 



,„) / ■ (2.7) 



where c*-^) is given by 



c(^) = , ' (2i 

^l + 4/OO^ 



and we have used (2.2) to eliminate Xj in the above. As a result, f^^> remains a function of 
(j)j whilst H, c^^\ V and 4>i generally depend on the collection of fields cj). Finally, we define 
the following slow variation parameters: 



H 



E^"' = E2»"'(f)^ 



f) = Y^^(ij) = Y^2c^''^^tLii 



Hi H ' 



J J 

The only general requirement on these parameters is e < 1, corresponding to a > 0, where a 
is the scale factor. We shall state explicitly when the additional restriction of slow variation 
is required, which corresponds to e^ ' ,7]^ ' ,s^ ' <^ 1. 

2.2 Perturbations and non-Gaussianity 

It is convenient to introduce the primordial curvature perturbation on uniform density hyper- 
surfaces C(t,x*), to characterise the scalar degree of freedom in the primordial perturbations 
(see [71, 72] for explicit definitions). The power spectrum P^ is then defined using the two- 
point function 

(Ck.Ck.) = {27rf Pdh) SH^i + ka), (2.10) 



where Ck is the Fourier transform of <^, k; are comoving wavevectors and 5^ is the three 
dimensional Dirac delta function. Moreover, the dimensionless power spectrum V(^{k) is 



often defined, related to the power spectrum by Vc_ = 2^P(;- Given this, the spectral index 
n^ characterises the deviation of V( from scale invariance 

dlogVc 

"C-1 = -JT^' 2.11 

a log K 

where scale invariance corresponds to n^ = 1. The two-point function completely defines 
the statistics of the field if C is purely Gaussian. Any signature of non-Gaussianity will be 
encoded in the connected contributions to higher order correlators, the next order being the 
three-point function 

(Ck.Ck.Cka) = {27TfB^{h,k2, h)5^iki + k2 + kg), (2.12) 



where the bispectrum B^^ is analogous to the power spectrum. We parameterise the deviation 
from Gaussianity by taking the ratio of the bispectrum to a combination of the power spectra 



6. ., , , ^_ B^{ki,k2,k3) 

pNL{k^,k2,ks) - p^^k,)P^^k2) + Pdh)Pdks) + Pdk2)Pdhy ^ ^^ 



where /nl is the A;-dependent non-linearity parameter^ and the numerical coefficient is a 
relic of the original definition in terms of the Bardeen potential $ = |^. Assuming a scale 

invariant dimensionless power spectrum and using the relation P^^ = ■^:^Pdk), the non- 
linearity parameter can be re-written as 



^/nl(A:i, A;2, k^) = ^-^ -^-^^ , (2.14) 



which is a more useful form for the following section. 

2.3 Super-horizon evolution and the 5N formalism 

The conversion of entropy to adiabatic modes during multiple-field inflation [15-17] allows 
the curvature perturbation Q to evolve on super-horizon scales [18, 19], and it is on this 
evolution we wish to focus. As such, we adopt the 6N formalism [18, 74-77], a powerful 
technique for evolving C, on super-horizon scales using only the field fluctuations at horizon 
exit and the homogeneous field evolution thereafter. 



''There are a number of sign conventions in the literature and we adopt those of [29] for ease of comparison. 
See [73] for a summary of the various conventions. 



In order to employ the 5N formalism we make two restrictions on the background 
dynamics. Since our scenario admits multiple sound horizons we demand that these are 
comparable whilst observable scales exit during inflation, such that c^^' = c^ for all / during 
this interval.^ Horizon exit^ therefore equates to evaluating a quantity when c^/c = flj^i/*. Our 
second restriction is that of slow variation at horizon exit, since this simplifies the spectrum 
of field fluctuations produced at horizon exit. 

Given these restrictions, the SN formalism uses the separate Universe approach [18, 75, 
76, 78] to allow the identification of the curvature perturbation (^ with the difference in the 
number of e-folds between the perturbed (A^) and homogeneous background (Nq) universes, 
evaluated between an initially flat hypersurface t^ (e.g. shortly after horizon exit) and a final 
uniform density hypersurface tf (e.g. early in the radiation dominated epoch) 



C{tf,x') = 6N{tf,x') = N{U,tf,x')-No{U,tf). (2.15) 



Assuming A^ is a local function of the scalar fields at horizon exit, which are split into a 
homogeneous background and local perturbation (j)^{ti,,x'^) = (/>^(i*) + (5^^(i*,x*), C can be 
expanded in powers of 6(/) as 

C = Y^ Nj6(^' + ^Y1 ^'IJ i^^'^^' - ^^^'^<t>')) + ■ ■ ■ , (2.16) 

/ IJ 



where Nj is with respect to the field / at horizon exit. Using this expansion, we are able 
calculate the relevant quantities (e.g. /nl) at time tf given the field fluctuations at time i^ 
and the homogeneous field evolution between t* and tf. For example, on substitution of the 
expansion (2.16) into the two-point function (2.10), one finds that the dimensionless power 
spectrum can be expressed as [77, 79] 



'Pc = yN^V,. (2.17) 



I 



Here we have defined the dimensionless power spectrum of scalar field fiuctuations at horizon 
exit using the two point function 

{HLM.) = (2vr)3 5^^|^P,<53(ki +k2), V. = {^^\ (2.18) 



^See [59, 60] for a discussion of when this is not the case. 

®For brevity we will refer to the 'sound-horizon' simply as the 'horizon'. There should not be any ambiguity 
since, in this context, the sound-horizon will be the only relevant scale. 



where we have used slow variation at horizon exit and 6^"^ is the kronecker delta symbol. 
The spectral index can be treated similarly, giving 

nc-l--2e^ + — ^^ — -2 . (2.19) 



Progressing to the three-point function, substitution of the expansion (2.16) into (2.12) neatly 
highlights the physical sources of non-Gaussianity 



UK 



I Y^ NiNjNKL{S<j)U^i^ {8<P'' * <5(/.^)k3) + 2 perms ) , (2.20) 



2 

IJKL 



where here • denotes a convolution and 'perms' denotes cyclic permutations over the mo- 
menta. The expression (2.20) shows two distinct contributions to the three-point function. 
The first term contributes to non-Gaussianity in (^ through the intrinsic non-Gaussianity of 
dip , produced by quantum field interactions on sub-horizon scales. Whilst this term van- 
ishes identically for Gaussian 6(p^ , C can still be non-Gaussian through the second term. This 
contribution is due to non-linear behaviour in (^ on super-horizon scales and is often referred 
to as local type.^ It is on this contribution that we shall predominantly focus. 

Neglecting the connected part of the four-point function and using Wick's theorem to rewrite 
the four-point functions as products of two-point functions, the latter term can be expressed 

as [77, 79] 

I Y, NjNjNKL{Scl>i,dct>i^{S^''*6ct>'^)^,) + 2perms = 



2 

IJKL 



(27r)3 Att'V^ ^i|i Y ^,iN,J^JJ '^'(ki + k2 + kg), (2.21) 



a^f ,. 



such that the bispectrum becomes 



.4^2E^ / 6 (3) ZijNjNjNjj 



B^{h,k,,h) = ^^'Vl^ -MikiMM + ^\' ''J,r ' (2-22) 



^See [80] for a discussion of the shape dependence of the bispectrum. 



and we have used (2.17) to replace "Pj, with "P^. Adopting the notation of [21], the /c-dependent 

(3) 

parameter f^-l{ki,k2,ks) accounts for the sub- horizon contribution for which we omit an 
exphcit expression, since we are primarily concerned with the super-horizon contribution. 

(3) 

We will return to /j^^ i^ our conclusions. Finally, on substitution into (2.14), we arrive at 

hL{ki,k2, ks) = /iS(fci, A;2, ks) + fl,^l (2.23) 



Again adopting the notation of [21], the A;- independent parameter^ f^-^ accounts for the 
super-horizon contribution and is given by 

eA^)_ EijNiNjNij , . 

5^NL- (^^Ar2^)2 ■ ^^-^4) 



Our aim then is to calculate P^, n,^ and f^^ in the model described above using the expressions 
(2.17), (2.19) and (2.24) respectively. To evaluate these expressions we must evaluate the 
terms Nj and Njj. Unlike single field inflation however, there is no unique attractor in the 
multiple-field case [83]. There are therefore an infinite number of trajectories and, as such, 
the derivatives of N become non-trivial. One is usually limited to numerical studies without 
a further restriction and so, to make analytical progress, we focus on scenarios satisfying a 
sum separable Hubble parameter in the following section. 



3 The role of rapidly varying sound speeds 

To investigate how local-type non-Gaussianity is affected by rapidly varying sound speeds, 
we apply the preceding results to a subset of analytically tractable cases. We note again 
that such behaviour was explored in the single field scenario by [64-68], in the context of 
equilateral type non-Gaussianity. In this section however, we show that rapidly varying sound 
speeds in the multiple-field case can, in principle, lead to large local-type non-Gaussianity. 

Given the lack of a unique attractor in multiple- field scenarios, certain restrictions are 
required to evaluate the expressions in section 2.3 analytically. For example, [21] calculated 
the non-linearity parameter analytically for canonical two-field models by demanding a sum 
separable potential y((/>, x) = V^^'{(l)) + V^^'{x)- This invoked the slow roll approximation 
throughout however, to reduce the second order field equations to first order expressions. 
To study more diverse dynamics, [29] developed an analogous formalism by writing the field 
equations in Hamilton- Jacobi form [43, 44], in which they are automatically first order, and 
considered a sum separable Hubble parameter H[(j),x) = H^'^'{(l)) + H^^'{x)- As a result, 
the slow roll approximation can be discarded after horizon exit and one can reliably study 
regimes in which slow roll is violated. 



See [81, 82] for scale dependent /nl- 



In this section we adopt the method invoked by [29] and demand a sum separable Hubble 
parameter. Whilst allowing the violation of slow variation after horizon exit, this also suits 
the case of non-standard kinetic terms, since the dynamics are more naturally described 
in the Hamilton-Jacobi formalism (see section 2.1). Given this restriction, we exploit the 
resultant integral of motion to derive analytic expressions for the derivatives of N and in 
turn V^, n^ and f^{^. The calculation closely follows [29] and we adopt analogous notation 
to clearly illustrate the additional effects of non-standard kinetic terms. 

3.1 Two-field model with a sum separable Hubble parameter 

We begin by considering the two-field case, with fields (j) and x, and demand a sum separable 
Hubble parameter 

H{4>,x) = H^^\^) + H^''Kx), (3.1) 



which leads to a number of simplifications. From equation (2.8), the sound speed c^ ' becomes 
a function of its respective field </>/ only, such that c^'^'{(j)) and c^^'{x)- Moreover, cross 
derivatives of H (i.e. H^^-^) become zero. The combined effect is to reduce the number of 
relevant slow variation parameters 



eW = 2cW('^y, 



,(x) = 2c(^) (^ 



2 



H H 

sW = 2c^^A, ,(x)=2cW^, (3.2) 



where e = e^'^^ -|- e^'^-' and we again emphasise that rj^^) and s^^' are not constrained to be 
small after horizon exit. The most crucial simplification however, is the ability to calculate 
the derivatives in section 2.3 analytically. To this end, we use (3.1) to write the number of 
e-folds between t^ and tj as 

N{U,tf)= Hdt = - w #+o / , , ^ dx, (3.3) 



which describes the evolution along the classical trajectory. Moreover, (3.1) permits the 
following dimensionless integral of motion 



+ / -^> (3-4) 



,9 "^ 

which is constant and unique for each classical trajectory, characterising motion orthogonal 
to it. By combining (3.3) and (3.4), we find the following expressions for the first derivatives 
of iV 






where, for clarity, we have included explicit references to tf and t* via the corresponding 
subscript. The task now is to calculate the derivatives of the fields at tf with respect to the 
fields at t^, where we can again invoke the integral of motion (3.4). 

3.2 Evaluating the field derivatives 

For brevity we focus only on -q^, with the remaining combinations following analogously. 
We begin by using the integral of motion (3.4) to write 

d^f _ dcpf dC 



dC 



(3.6) 



With respect to the first term, we note that since tf is defined as a constant density hyper- 
surface, H{(j)f,Xf) is a constant in a spatially flat Universe. As a result 



'iHf + ^h(x) = 0. 



dC '^f dC '^f 
Furthermore, differentiating (3.4) along C gives 



(3.7) 



^/ ■'^<^!'/ ■-/ ^ 'X/ 

which, combined with (3.7) yields the first term in (3.6) 

- 10- 



1 d(pf 1 dxf 

^^^(0) dC ^ c^x)jj{x) dC ^^' ^^' 



dC 



rrW f 1 , 1 



(3.9) 



The second term is simply given by the definition of C 

dC 1 



>)^W' 



(3.10) 



and we combine (3.9) and (3.10) to find the result for -^^ 

dcpf _ 1 



(3.11) 









Finally, we use the slow variation parameters (3.2) and include the remaining combinations 
of derivatives to find 






(3.12) 



which are used to evaluate the derivatives of A^ in the following section. 

3.3 Results 

Given the expressions (3.5) and the results (3.12), the first derivatives of A^ can be calculated 
in terms of slow variation parameters 



N^^ = ^ u, N^^ = ^ V. (3.13) 
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For brevity we have introduced the following definitions 



u= * \ v= * \ Zf = ^-^ ^^. (3.14) 



The above neatly demonstrates the nature of the 5N formalism, since the term Zf controls 
the super-horizon evolution whilst the remaining terms are evaluated at horizon exit. To 
calculate n^ and f^^ we also require the second derivatives of A^, which are found by simply 
differentiating (3.13) to give 






iVv.v, = - ^^ - -^ — . . . . ^ V ^'^* 



^^*^* 2 uw Ax)(x) ,/i^^r^ H, 



-* c* L.* -i/ei'-'c; 






where the derivatives of Zf are given by 



z,._= ^"' - i Z,y, = ^E^A. (3.16) 

0^ ,/;i^ 



Again for brevity we have introduced the following definition: 



2 ^(0)^(X) /i ^ss , ^ss 



tt2 AVI ,\XI / 1 „ss -I 

^/'/ ^/ /I ^/ 1 



•^-#^l2-t-2l7l- *^-^^> 



This is equivalent to the parameter A in [29] and both expressions contain the second term 
in parenthesis 

^ss ^ <^_Ji ^^ '/ , (3.18) 



12 



We find an additional term however, specific to non-standard kinetic terms, given by 



s'' = ' ^ +' ^ ■ (3.19) 



(2.19) and fi^l (2.24) to find 



Finally, we substitute the results for N j (3.13) and N jj (3.15) into the expressions for n^^ 

f(4) 
'NL 



ull- ^'* f^/ ' u ] +v[l- ^'* ,^,* ' v 



n^-l = _2e, - 2^ '- -^— ^ ^ ^, (3.20) 



u^ ^^ 






-J^l = — ^ T^. — • ^^-^^^ 



The above results are valid for the two-field scenario detailed in section 2.1 assuming com- 
parable sound speeds and slow variation at horizon exit, in addition to a separable Hubble 
parameter. Note that c; and c;'^ do not appear explicitly in these expressions, since we 
assume they are approximately equal and therefore cancel. 

It is straight forward to confirm that (3.20) and (3.21) recover previous results in the 
appropriate limits. To recover the single field expressions we consider x ~^ 0, such that 
Zf — )• Hf , u — >• 1 and f — >• 0, since in the single field limit the curvature perturbation C is 
conserved after horizon exit and all parameters can be evaluated at i*. With this, we find 



nc - 1 = -44^) + 2,i^) + 24^), lit = ei^) - ,i^^ - si^\ (3.22) 



and recover the single field results [45]. Moreover, it is trivial to check that we recover the 
canonical results, since in the limit c^'^'= o-^>= 1 the expressions (3.20) and (3.21) coincide 
with those of [29]. 

With regard to the general behaviour of our results, inspection of (3.20) shows that 
n^ — 1 will remain 0{e^), where here e represents first order in slow variation. This is in 
keeping with observation of a nearly scale invariant power spectrum [84]. Notice that whilst 
the above does contain additional factors of si- and s^^ with respect to [29], their effect 
does not alter the conclusion that uq — 1 remains 0(e*). 
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Turning to f^-^, the first two terms in parenthesis in (3.21) are similarly O(e^) and so 
would contribute little towards any enhanced non-Gaussianity. Moreover, barring certain 
fine-tuned choices for the trajectory, the prefactor of A is generally 0(1) (see [27] for a 
summary of when this term is large in the canonical case). As a result, A provides the only 
potential source of significant non-Gaussianity. Whilst the terms preceding the parenthesis 
in the expression for A (3.17) are 0(e*), tj^ and s^ can become much larger than unity. In 
the canonical case for example, [29] show that 7]^ can produce enhanced non-Gaussianity 
as the trajectory descends a potential in which r]^'!'' and rj^^' become large. For non-trivial 
sound speeds however, we find the additional factor s^r' which can, in principle, produce a 
similar contribution. For example, one might expect an analogous scenario to [29] in which 
■s"^** produces significant non-Gaussianity as the trajectory traverses a warped region in which 
c^'^' and c^^' change abruptly and so s^'^' and s^^' are enhanced. In the following section we 
will apply our general formulae to such a scenario and demonstrate how s^ becomes large, 
producing enhanced local non-Gaussianity during inflation. 

We note here that the expressions (3.20) and (3-21) are well suited for describing the 
evolution of n^ and f^^ during inflation. However, these are not necessarily the observed 
values since in general we must track their evolution from the end of inflation until they are 
imprinted on the cosmic microwave background (CMB) at decoupling. Ideally, one would 
know the entire evolution history between these times but, given our lack of knowledge of 
the early Universe, this is not currently feasible. In light of this, recent work has considered 
whether such non-Gaussianity produced during inflation can imprint upon the CMB [33-40] . 
For example, [33, 34] consider scenarios in which a limit is reached during or shortly after 
inflation, in which perturbations become essentially adiabatic and approximately conserved. 
It may be that the potential has a focusing region in which trajectories converge before the 
end of inflation and, in certain circumstances, can imprint upon the CMB. Alternatively, a 
waterfall field may end inflation whilst f^^ is large, after which it is preserved. If a conserved 
limit is reached shortly after inflation however, one must still choose a method for reheating 
and the results may depend upon this choice. As a further example, [35, 36] consider sum 
separable potentials and find non-Gaussianity is damped during the transition to single field 
behaviour, leaving little imprint on the CMB. Such study is invaluable for understanding 
what we can infer about inflationary dynamics from observations of non-Gaussianity and it 
would be interesting to make similar considerations in our case. In this work however, we 
simply illustrate the potential production of non-Gaussianity through multiple-field dynamics 
during inflation, facilitated by non-trivial sound speeds as opposed to the potential. 



4 Inflation in two cutoff throats 

To illustrate the more general arguments presented in the previous section, we present a 
specific scenario that can be described by the above framework. In particular, we consider 
two-field inflation in the tip regions of two warped throats. We find that the combined effect 
of rapidly varying sound speeds and multiple-field dynamics can lead to significant local 
non-Gaussianity towards the end of inflation, for certain initial conditions. 
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Figure 1. Geometric illustration of our example model. Inflation is driven by two probe D3 branes 
traversing two distinct warped throats glued to a compact Calabi-Yau manifold in type IIB string 
theory. The radial co-ordinates of the D-branes play the role of the fields 4> and %. The throats are 
produced by stacks of D3 branes that source RR fluxes and we include fractional D3 branes at the 
tips to cut-off the throats. 



4.1 Background trajectory 

We model inflation as driven by two probe D3 branes traversing two distinct warped throats 
glued to a compact Calabi-Yau in type IIB string theory [58-60], as in figure 1. Each throat 
is produced by a stack of N D3 branes located at the tip of the conifold, sourcing N units of 
5-form flux which warp the space. In addition, we add M wrapped D5 branes at the apex 
of each throat that source M units of magnetic RR 3-form flux, through the 5^ of the basis 
T^^ of the cone. The resulting space-time is a Klebanov-Tseytlin throat [85] and we refer to 
[86] for a review of the geometrical aspects of this set up and [87] for the case of single field 
inflation in a cut-off throat. Given the above, the warp factors are given by 



''"-j'i^-'^Mi 



(4.1) 



where f^^' is given by replacing (/> — )■ %. For simplicity we have assumed the same warping in 
each throat, such that A^ = A^'^ = Aj, which does not qualitatively alter our conclusions. 
Here Aj control the number of RR fluxes switched on and, in particular, A2 is dependent on 
M [86] . Geometrically, the logarithmic running of the warp factor implies that the effective 
number of units of 5-form flux varies, depending on the position along the throat. Note 
the infrared singularity at (/> = Ase"^''^^ where the supergravity approximation used to 
obtain the geometry breaks down and the solution becomes unreliable. By considering a 
generalization of the previous set-up, obtained by wrapping branes on a deformed conifold, 
this singularity can be tamed [88]. However, the inflationary trajectory we consider does 
not reach the IR singularity and so we adopt (4.1) for the warp factors. Figure 2 illustrates 
jy<p) for representative values of the parameters, which we will use hereafter. Note that such 
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Figure 2. Left: The warp factor /^"^^ (solid) and corresponding sound speed c^'^^' (dashed) plotted 
as functions of (/), with Ai = 6 x 10l^ A2 = 2, A3 = 1 and U'-'^^ = 1.188 x 10"^. Right: The rapid 
increase in c^*^) towards the tip of the throat in further detail. 



choices are used for illustrative purposes only, since a full exploration of parameter space is 
beyond the scope of this work. Given this, we choose a linear, separable Hubble parameter 



H{(I),x)='H'-^U + 'H'-''^ 



X, 



(4.2) 



such that rj^'*''= r]^^>= 0. This guarantees that the only significant contribution towards f^-^ 
in (3.21) must be s**, since ry^* = throughout inflation. Finally, with expressions for the 
Hubble parameter and warp factors, the potential is fixed by the Friedmann equation (2.7) 
and, throughout most of inflation, is well approximated by 



Vi(f>, x) ^ 3?^^^^ + 3?^W'x^ + 6?^W^(>^Vx- 



(4.3) 



Whilst this potential is required to suit our analytical approach, it allows us to demonstrate 
properties that we expect to be representative of a broader class of potentials. Using (2.8), 
we arrive at the following expression for the sound speeds 



M 



l + ^^^^(l + A2log(^ 



A3 



(4.4) 



where c^^' is again given by replacing — t- %. Figure 2 illustrates c^"^' for the f^'^' plotted pre- 
viously. This provides some qualitative insight into the inflationary dynamics. Considering 
for a moment the single field case, infiation progresses efficiently with c^'''' <C 1 and e*-*^^ <C 1 
for (j) > 0.7, as in the standard single field DBI case. For (j) < 0.7 however, corresponding 
to the tip of the throat, the sound speed rises rapidly and drives e^'''' towards unity, where 
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Figure 3. Left: The trajectory in field space, originating at 0(t*) = x(^*) — 1 ^-^id ending after 
A'' ~ 60 e- folds of inflation, when e = 1. Right: Enlarged region of the trajectory (solid) illustrating 
the turn in the x direction towards the end of inflation. The dashed line shows the straight line 
trajectory corresponding to "H^"^' = T-L^^i for comparison. 



inflation ends. Similar qualitative arguments hold for the two-field case but the evolution 
will depend more on the choice of model parameters. 

To study the precise dynamics we must solve the equations of motion. Since the Fried- 
mann equation is automatically satisfied we need only solve the field equations (2.6) for a 
given choice of parameters and initial conditions. As such we choose Ai = 6 x 10^^, which is 
typically required in standard DBI [49], A2 = 2 and A3 = 1 to characterise the warp factors. 
Furthermore, we choose (/>(t*) = x(^*) = 1 as our initial conditions, where t^, is the time 
at which observable scales exit the horizon. We require A^ ~ 60 e-folds between t^ and the 
end of inflation, which we choose to define as e = 1. Note that this choice is a working 
definition for the end of inflation, since we do not consider an explicit model of reheating 
in this work. The choice of parameters and initial conditions thus far is 'symmetric', in the 
sense that for T-L^'^' = T-L^^' the resultant trajectory would be a straight line. To facilitate 
multiple-field behaviour we introduce a small asymmetry such that T-L^'^' = 1.188 x 10"^ and 
y[x) = 1.192 X 10~^. With this choice of parameters c^'f'^ /c^^^ ~ 1 + 10"^ at U and is therefore 
consistent with our approximation that the sound horizons are comparable when observable 
scales exit. Given this, we solve the field equations and plot the infiationary trajectory in 
figure 3. 

As we would expect, given the almost symmetric initial conditions and parameter 
choices, the trajectory is approximately straight. However, towards the end of infiation 
the trajectory curves in the x direction. To understand this, we show the evolution of the 



sound speeds parametrically in figure 4. Throughout most of infiation c' 



.(0) 



c(^) <C 1, be- 



having essentially as standard single field DEL At the tip of the throats, corresponding to 
the last few e-folds of infiation, both sound speeds increase rapidly. By choosing 71^^' > T-L^'^' 
however, c-^' increases slightly before c-'^', producing the upwards turn in figure 4. Since the 
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Figure 4. Evolution of tlie sound speeds c'"^) and c^^'' (solid). For most of inflation 



U) 



J.X) ^ 1 



but during the last few e-folds both rise rapidly, with c^^i > c^'^'. This explains the downward 
curving trajeetory in figure 3. For eomparison, the dashed line shows the corresponding evolution for 

7^(0) ^u^x)^ 



motion of x is now less inhibited than that of (p, the trajectory in figure 3 curves in the x 
direction. Whilst this implies the presence of isocurvature modes at the end of inflation, we 
cannot track any subsequent evolution without an explicit model of reheating, as noted in 
section 3.3. 



4.2 Evolution of /, 



(4) 

NL 



Given this trajectory, we use the expressions in section 3.3 to study the evolution of the 
relevant quantities (e.g. /^l) as a function of tj for a fixed t^. Considering first the two- 
point statistics we find V(^ = 2.44 x 10~^ and n^ — 1 = —0.0108 at the end of inflation, 
which are consistent with observations [84] and remain approximately constant throughout 
inflation. This is as expected, since the corresponding expressions are not sensitive to the 
rapidly varying sound speeds at the end of inflation. 

Considering now the three-point statistics, we find that the rapidly changing sound 
speeds allow s^'^' and s^^' to increase towards the end of inflation. We note again that, since 
our approach does not rely on slow variation after horizon exit, we are able to study such 
regimes with confidence in our results. We therefore find that s** increases dramatically, 
as illustrated in figure 5, and so therefore does A. In addition to the enhancement of A, 
the small departure from single field behaviour towards the end of inflation, as illustrated in 
figure 3, produces a non-zero term preceding A in the expression for /j^^ (3.21). Note that 
regardless of the magnitude of A, /j^^ would remain small for a straight trajectory (i.e. for 



-20 at the end 



j{{<t>) = 'H(^)). The combined effect therefore is a rapid increase towards f^^ ~ 

of inflation, as shown in figure 5. Note that the negative sign of the non-linearity parameter 
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Figure 5. Left: Evolution of s'^'^ with respect to tf, plotted as a function of e for the trajectory in 
figure 3. The rapidly varying sound speeds produce s^* ~ —34000 at the end of inflation, when e — 1. 
Right: The corresponding evolution of fj^-^ , where the combination of rapidly varying sound speeds 
and a turn in the trajectory produces /j^^ — ^20 at the end of inflation. 



in this example is the result of rapidly increasing sound speeds, which is disfavoured but not 
excluded by observations [84]. It would be interesting to apply similar techniques to [33, 34] 
to further understand the sign of /j^^ ™ such scenarios. 

It is important to note that DBI type models are usually considered in the context of 
equilateral type non-Gaussianity, through enhanced field interactions at horizon exit. For 



example, for single field DBI in the equilateral configuration one finds |/j 

(3) 



(3), 

nlI 



c„ ^ at horizon 



exit, which is preserved thereafter [49]. Moreover, /nl — f^l to a very good approximation, 
since the local contribution is negligible in the single field case. This places a firm lower 
bound on the sound speed at horizon exit. In the scenario considered above the sound speeds 
at horizon exit are similarly small and the naive expectation would be for a prohibitively large 
/j^L, in addition to the large /j^^ '^^ have already discussed. The introduction of multiple field 
dynamics can alter this conclusion however, via conversion between entropy and adiabatic 
modes both during horizon exit and thereafter (see, for example, [50-52, 61] for the case 

(3) 

of multiple field DBI in the context of a single throat, in which f^^ is suppressed by such 
dynamics). We leave an explicit calculation of the equilateral contribution and its correlation 
with the local contribution to a follow-up paper and, for now, emphasise our main result of 
large local type non-Gaussianity from dynamics associated with non-trivial sound speeds. 

Finally we note that, given our working definition of e = 1 as the end of inflation, at 
which time we evaluate the relevant quantities, we are unable to track their evolution from 
the end of inflation until they imprint upon the CMB. Whilst the implementation of an 
explicit model for reheating or the imposition of a transition to adiabaticity is outside the 
scope of this work, it would be an interesting topic for future research. 
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Figure 6. Left: Evolution of ^nl with respect to tf, plotted as a function of e for the trajectory in 
figure 3. Right: log I^nlI as a function of log |/nl|- The result is a straight line with slope ~ 2, such 
that gNL oc /^L i^ this case. 
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Figure 7. Left: Evolution of tnl with respect to tf, plotted as a function of e for the trajectory in 
figure 3. Right: log |tnl| as a function of log |/nl| (solid). The result is a straight line with slope 
~ 1.7, such that tnl oc f^'^ in this case. The Suyama-Yamaguchi inequality (dashed) is also shown. 



4.3 The trispectrum 

In addition to the bispectruni, it is interesting to briefly consider the behaviour of the Fourier 
transform of the four-point function, the trispectrum. In particular, we are interested in 
the local- type ^nl and tnli the analogues of f^{^ for the trispectrum. By performing an 
analogous calculation to that in section 2.3, (^nl and tnl can be recast using the 6N formalism 
[79, 89, 90] 



9NL 



25 EijK NjjkNjNjN^k 



54 



(EiNf) 



2\3 



TNL 



Zijk^JjNikNjNk 



(4.5) 
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It is simple to evaluate tnl in principle, since we already have expressions for the first and 
second derivatives of A^. Furthermore, we need only take a further derivative of (3.15) to find 
the required terms in (7nl- In practice however, such expressions are prohibitively long. As 
such, we restrict ourselves to presenting their evolution along the inflationary trajectory. In 
figure 6, (7nl oc f^^ and becomes O(IO^) towards the end of infiation. Similarly, tnl becomes 
0(10^) in figure 7, where tnl oc f^'^. Note that the Suyama-Yamaguchi inequality [91, 92], 

7"NL ^ ( f /nl ) ' remains satisfied by a large constant of proportionality. 



5 Conclusions 

We have studied the effect of non-trivial sound speeds and multiple-field dynamics on non- 
Gaussianity during inflation. In particular, we have shown that rapid changes in the sound 
speeds can, in principle, produce local-type non-Gaussianity during a turn in the trajectory. 

As an example of multi-component inflation with non-trivial sound speeds, we consid- 
ered a multiple-DBI model and used the 6N formalism to track the super-horizon evolution of 
perturbations. Having expressed the homogeneous equations of motion in Hamilton-Jacobi 
form, we adopted the method of [29] and used a sum separable Hubble parameter to de- 
rive analytic expressions for the relevant quantities in the two-field case, valid beyond slow 
variation. 

We find that non-trivial sound speeds can, in principle, produce significant local-type 
non-Gaussianity. Deviations from slow variation, such as rapidly varying sound speeds, 
enhance this effect. Note that this contribution is distinct from that of the potential and 
the general result is a combination of both sources. To illustrate our results we considered 
infiation in the tip regions of two warped throats, in which the sound speeds rapidly increase 
at the end of inflation. This, alongside a turn in the trajectory, produces large local-type 
non-Gaussianity towards the end of inflation. The non-linearity parameter is negative in this 
example however, owing to the rapidly increasing sound speeds. 

There are a number of outstanding questions that lend themselves to further work. 
Most noteworthy is that, in this paper, we considered the local-type contribution to /nl 
from multiple-fleld dynamics during inflation. It would therefore be interesting to compare 
this to the equilateral contribution produced at horizon exit by small sound speeds. Whilst 
this contribution is dominant in the single held case, the introduction of multiple fields may 
alter this conclusion via conversion between entropy and adiabatic modes (see, for example, 
[50-52]). Moreover, correlations between different contributions may lead to one of the few 
explicit examples of models characterised by combined equilateral and local non-Gaussianities 
(see [61] for an alternative). We intend to address these issues in much greater detail in a 
follow-up paper. 

Other open questions remain regarding the application to more general cases. Firstly, 
whilst our example trajectory illustrates the production of local-type non-Gaussianity through 
non-trivial sound speeds, it is unclear how generic this behaviour is. It would therefore prove 
useful to study a wider parameter space (see [27] for an example in the canonical case). Sim- 
ilarly, our use of a sum separable Hubble parameter, whilst providing analytically tractable 
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expressions, also restricts the choice of potential. Thus it would be interesting to study 
more general cases with alternative analytical or numerical methods. Note that this is also 
true in the canonical case. Finally, we highlight that our model based approach explicitly 
demonstrates the dynamics necessary for local-type non-Gaussianity. This must be consid- 
ered in light of more general methodologies however, such as an effective field theory approach 
[93, 94], which offer wider applicability. It would therefore be interesting to bridge the gap 
between the two. 

We end by again noting that our approach treats multiple-field dynamics during inflation 
and, since we did not consider a model of reheating or an approach to adiabaticity, we 
cannot conclude that such values are necessarily observed in the CMB. In our illustrative 
example however, the non-linearity parameter peaks at the end of inflation, corresponding to 
a violation of slow variation at the tips of the throats. We would argue that such scenarios 
seem the most likely to produce non-Gaussianity during inflation that can imprint upon the 
CMB. More generally however, it would be interesting to apply similar techniques to [33-40] 
in the case of non-standard kinetic terms, to further understand what we can infer about the 
dynamics of inflation from observations of non-Gaussianity. 
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A Equations of motion in the Hamilton-Jacobi formalism 

In this appendix we derive the equations of motion (2.2) and (2.3), presented in section 2.1, 
in Hamilton-Jacobi form. This essentially extends the calculation in [43] to non-canonical 
cases satisfying the following action 



iyd*i-V^[R + 2^P,-2F], (A.l) 



*=2 



where R is the Ricci scalar and g is the determinant of the metric tensor g'^^^ . Note that Pj 
is a function of the single scalar field (pj and kinetic function Xj = —29'^'^4>i,fj.4>i,u, whereas 
the potential F is a function of the set of scalar fields cj) = {<j)i,(j)2, ■■■,4'n}- Note also that 
we return to the canonical case with Pj = Xj. 

It is useful to recast the action (A.l) using the Arnowitt-Deser-Misner (ADM) formalism 
[69], originally conjectured as a Hamiltonian formulation of general relativity. The formalism 
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proves useful for studying the background evolution as it naturally produces scalar field 
equations in Hamilton- Jacobi form [43, 44], which is more suited to non-trivial sound speeds. 
Moreover, the formulation introduces the lapse function N and shift vector N^ which under 
variation become Lagrange multipliers, simplifying the equations of motion for perturbations. 
The ADM metric is given by 



ds^ = -N^dt^ + hij {dx' + N'dt) {dx^ + N^dt) , (A.2) 



where hij is the spatial 3-metric. In terms of this metric, the action (A.l) and kinetic term 
becomes 

S =- I d'^xVh[NR^'-^'> + NKijK'^ - NK^ + 2iV ^ P/ - 2NV] , (A.3) 

^/ = ^ {<Pi - A^V/,i)' - l^Pi^i'P'i, (A.4) 



where K = K^, R^^' is the three dimensional Ricci scalar and indices are raised and lowered 
using the spatial metric. Kij is the extrinsic curvature, given by 



where u denotes the covariant derivative with respect to the spatial metric. It is convenient 
to define the traceless part of a tensor with an overbar and we use the extrinsic curvature as 
an example 

Kij = Kij - h,jK. (A.6) 



It is trivial to check that Kij is indeed traceless. Variation of the action (A.3) with respect 
to A^ and A^* yields the energy and momentum constraints 

KijK'^ - '^K^ - ii(3) + 2p = 0, (A.7) 

^i-^^K + En^'^Ai = o. (A.8) 

Notice that these relations are identical in form to the canonical results of [43] , since the effect 
of the non-canonical kinetic terms are packaged into the energy density p and momenta 11/ 
of the fields 
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n/ = ^ (0/ - ivV/,0 ■ (A.io) 



Similarly, variation of (A. 3) with respect to the spatial metric yields the gravitational field 
equations 

K - N^K\, + ivi; - N {^^K,jK'^ + 1^2 + 1^(3) ^ 1 5^ ^ 0^ (^.11) 

1 -H^A* _ AT fl^i^i ^ r(3)« 



-iV|';5] - AT (^KKj + R'"''^ - S}) = 0, (A.12) 

where the stress three tensor is defined as 



Sij = Y, ( Pi,Xj (^I,^^I,j + hij Pi j - hij V. (A. 13) 



Again the non-canonical contributions are packaged into (A. 13). Finally, the field equations 
follow from variation with respect to (pi 

NPij -NVi-^ ( V^P/,x,n,) + (p,,x,iV'n,) + (^NPi,Xj(l)'i) = 0. (A.14) 



In the canonical limit the above system of equations recover the results of [43]. 

To derive a set of equations for the homogeneous evolution in Hamilton-Jacobi form we 
again follow [43] by expanding in spatial gradients. We consider only terms up to first order 
in spatial gradients (e.g. neglecting terms such as A'', and R ) and set A^* = in order to 
simplify the system of equations. The analysis follows [43] very closely and so for brevity we 
simply quote the results. The momentum constraint (A. 8) becomes 

H,, = -^Y.Ui^i,„ (A. 15) 

I 



where we have replaced the trace of the extrinsic curvature with the spatially dependent 
Hubble parameter H[t,x'^) = g^ . Since Hi = "^i Hjcpi^i, matching terms with (A. 15) 
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and considering a spatially flat Friedniann-Robertson- Walker Universe (such that A^ = 1) 
gives the field equations 

Hj. (A.16) 



Pi,x, 



The Friedmann equation follows directly from either the energy constraint (A. 7) or the equa- 
tion for the trace of the extrinsic curvature (A. 11) 



/ 4jj2 \ 



Finally then, we arrive at the equations of motion (2.2) and (2.3) used in the beginning of 
section 2.1. 
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